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ABSTRACT 


This  study  was  conducted  to  examine  the  effects  of  clouds  on  the  ocean 
mixed  layer,  both  short-term  and  seasonal.  It  utilized  data  collected  at  Ocean 
Station  Papa  in  the  northeast  Pacific.  Two  numerical  modeling  simulations 
were  performed  (i.e.  with  variable  cloud  and  with  variable  precipitation). 
The  results  for  the  variable  cloud  simulation  indicated  that  the  downward 
surface  buoyancy  flux  and  longer  daylight  period  in  summer  may  induce  a 
significant  albedo  effect  of  cloud  on  ocean  mixed  layer.  The  upward  surface 
buoyancy  flux  and  longer  night  period  in  winter  will  result  in  a  pronounced 
greenhouse  effect  of  cloud  on  ocean  mixed  layer.  The  results  of  variable 
precipitation  simulation  showed  that  the  mixed  layer  is  most  sensitive  to 
precipitation  between  October  and  March. 

Model  predictions  are  verified  using  data  at  Ocean  Station  Papa  for 
monthly  and  yearly  mean  values  of  cloud  cover  and  precipitation.  The 
comparison  between  model  prediction  and  observations  shows  that  the  mean 
values  of  observed  MLD  (  H  =  60.9  m)  are  much  deeper  than  model-predicted 
values  (  H  =  36.5  m). 
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INTRODUCTION 


In  general,  both  clouds  and  precipitation  can  have  a  significant  influence 
on  upper  ocean  dynamics.  Alteration  of  short  wave  and  long  wave  radiation 
at  the  ocean  surface  caused  by  clouds  and  variation  of  upper  ocean  salinity 
caused  by  precipitation  and  evaporation  may  also  change  the  ocean  mixed 
layer  structure.  Very  few  one-dimensional  mixed  layer  experiments  had 
previously  explained  how  the  clouds  and  precipitation  affect  the  ocean  mixed 
layer.  It  might  be  that  the  cloud  amount  and  precipitation  were  inaccurately 
observed  or  measured,  or  that  the  water  mass  flux  has  been  given  less 
attention  by  oceanographers,  although  it  is  realized  that  clouds  and 
precipitation  contribute  to  the  surface  buoyancy  flux  that  influences  depth 
and  temperature  in  the  ocean  mixed  layer. 

It  is  the  intention  of  this  study  to  examine  the  effects  of  clouds  on  the 
evolution  of  the  gross  thermodynamic  properties  of  the  ocean  mixed  layer. 
Employing  data  derived  from  observations  taken  at  Ocean  Station  Papa  (50°N, 
145°W)  we  shall  focus  specifically  on  the  short  and  long-term  fluctuating 
relationship  between  degree  of  cloud  cover  and  its  oceanic  consequences. 

Earlier  studies  of  ocean  mixed  layer  dynamics  dealt  explicitly  with 
equations  for  the  production,  alteration  and  destruction  of  turbulent  kinetic 
energy  within  the  mixed  layer.  Kraus  and  Turner  (1967)  were  the  first  to 
examine  the  turbulent  kinetic  energy  budget  in  a  one-dimensional  mixed 
layer  model  for  the  ocean,  using  the  approximately  decoupled  equations  for 
the  thermal  and  mechanical  energy.  It  was  assumed  that  the  heat  input  at  the 
ocean  surface  was  known,  and  the  water  mass  entrainment  at  the  base  of  the 
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mixed  layer  needed  to  be  predicted.  The  wind  stress  is  treated  as  a  production 
term  that  generates  the  turbulent  kinetic  energy  for  mixing.  The  entrainment 
was  associated  with  wind  energy  input  to  the  water  column  minus  the  work 
performed  in  overcoming  the  buoyancy  flux  throughout  the  mixed  layer. 
The  problem  of  the  Kraus-Turner  model  was  that  it  only  considered  the  heat 
flux  and  wind  forcing  and  neglected  the  viscous  dissipation  and  salinity 
effects.  The  dissipation  was  assumed  to  be  a  fixed  fraction  of  wind  stress 
production  in  the  model  of  Geisler  and  Kraus  (1969),  Miropol'skiy  (1970)  and 
Denman  (1973).  The  further  parameterization  of  dissipation  is  needed  in 
certain  instances  (Elsberry  et  al.  (1976),  and  Kim  (1976)). 

The  first  study  which  considered  salinity  and  its  effects  on  the  density 
structure  in  the  ocean  mixed  layer  was  by  Miller  (1976).  In  his  study,  he 
indicated  that  without  salinity  effect,  the  cooling  induced  by  entrainment  was 
greater  than  the  heating  from  the  surface,  resulting  in  ocean  mixed  layer 
cooling.  On  the  other  hand,  with  salinity  included  the  layer  deepens  more 
slowly  and  entrains  less  cold  water  at  the  bottom  allowing  surface  heating  to 
dominate  which  causes  the  mixed  layer  temperature  to  increase. 

A  one-dimensional  bulk  model  of  ocean  mixed  layer  which  included 
dynamic  instability  with  turbulent  erosion  was  introduced  by  Garwood  (1977). 
This  model  modulates  the  mean  entrainment  rate  by  diurnal  heating  and 
limits  maximum  layer  depth  by  enhancing  the  effect  of  viscous  dissipation. 
His  study  also  included  the  possibly  important  effect  of  long-term  salinity 
fluxes  on  the  mixed  layer  thermal  profile  due  to  surface  heat  fluxes, 
precipitation  and  entrainment  at  the  base  of  the  mixed  layer.  This  model  was 
chosen  for  this  study. 
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Livezey(1988)  used  an  oceanic  mixed  layer  model  to  examine  the  effects  of 
discrete  precipitation  events  on  the  short-term  and  seasonal  evolution  of 
ocean  mixed  layer  temperature  and  salinity  structure  at  Ocean  Station  Papa. 
He  indicated  that  the  effect  of  precipitation  on  ocean  mixed  layer  largely 
depends  on  the  strength  of  wind  forcing.  When  wind  speed  is  low,  the 
fluctuation  of  mixed  layer  temperature  and  salinity  will  be  considerably 
enhanced  by  the  precipitatioin  The  model  results  suggest  that  the  amount  of 
precipitation  observed  at  Ocean  Station  Papa  was  too  low  to  explain  the 
observed  ocean  salinity  structure  for  the  year. 

The  importance  of  cloud  cover  on  the  mid-latitude  ocean  mixed  layer 
has  not  been  investigated  in  a  numerical  model,  although  it  plays  an 
important  role  in  heat  budget  of  the  ocean  surface.  Recent  work  about  the 
thermodynamical  process  between  cloud  and  tropical  ocean  mixed  layer  was 
introduced  by  Chu  and  Garwood  (1989,  1990).  They  show  that  the  ocean 
mixed  layer  and  clouds  are  coupled  by  both  heat  and  moisture  fluxes  across 
the  air-ocean  interface.  In  this  coupled  system,  as  far  as  the  turbulent  kinetic 
energy  is  concerned,  the  fresh  water  influx  at  the  ocean  surface  due  to  excess 
precipitation  over  evaporation  is  a  damping  factor  and  the  surface  cooling 
caused  by  reduction  of  solar  radiation  at  the  ocean  surface  by  cloud  is  a  forcing 
factor.  The  aim  of  this  study  is  to  examine  the  effects  of  these  two  factors  on 
the  mid-latitude  ocean  mixed  layer  dynamic  by  using  meterological  and 
oceanographic  observations  in  1959  at  Ocean  Station  Papa. 
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n.  EFFECT  OF  CLOUD  ON  OCEAN  MIXED  LAYER 


Generally  speaking,  there  exist  two  thermodynamic  processes  in  the 
coupled  cloud-ocean  mixed  layer  system.  First,  the  incoming  solar  radiation 
at  the  ocean  surface  will  be  reduced  by  scattering  and  absorption  of  cloud, 
which  cools  the  ocean  surface  layer  and  increases  mixed  layer  entrainment. 
This  is  a  forcing  factor.  Second,  if  rainfall  begins,  the  precipitation  will  dilute 
the  surface  salinity,  stabilizing  the  upper  layer  and  decreasi.ig  mixed  layer 
entrainment.  This  is  a  damping  factor.  Figure  1  shows  the  relation  between 
the  clouds  and  ocean  mixed  layer  (Chu  and  Garwood  1989). 


Figure  1.  The  Relationship  between  Clouds  and  Ocean  Mixed  Layer 
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It  should  be  mentioned  that  deepening  or  shallowing  of  the  ocean  mixed 
layer  depends  on  the  ocean  surface  buoyancy  flux  and  surface  wind  stress. 
Wind  stress  always  generates  turbulence,  tending  to  deepen  the  mixed  layer. 
However,  the  buoyancy  flux  has  two  different  effects.  Upward  buoyancy  flux 
make  the  v’'ater  column  statically  unstable,  which  generates  turbulence  and 
deepens  the  ocean  mixed  layer.  Downward  buoyancy  flux  makes  the  water 
column  statically  stable,  which  damps  the  turbulence  and  possibly  shallows 
the  ocean  mixed  layer.  When  the  ocean  mixed  layer  deepens,  the  mixed  layer 
entrains  water  from  below  and  increases  the  potential  energy.  If  there  is 
relatively  stronger  downward  buoyancy  flux,  the  shear  production  may  be 
insufficient  to  overcome  the  buoyant  damping  and  turbulence  cannot  reach 
the  former  mixed  layer  depth,  causing  mixed  layer  retreat. 

In  the  diurnal  period,  the  clouds  also  have  two  important  effects  on  the 
radiation  balance  at  the  ocean  surface:  namely,  an  albedo  effect  when  the 
downwelling  solar  radiation  is  decreased  due  to  reflection  and  absorption  of 
incoming  solar  radiation  by  clouds  during  the  daytime,  and  a  greenhouse 
effect  when  the  upwelling  longwave  radiation  is  decreased  due  to  the 
emission  of  long  wave  radiation  back  to  the  ocean's  surface  by  clouds.  The 
albedo  effect  on  ocean  mixed  layer  is  dependent  on  diurnal  and  seasonal 
variations  of  solar  radiation.  During  summer,  solar  radiation  is  stronger, 
therefore  the  albedo  effect  of  cloud  significantly  impact  the  ocean  mixed  layer. 
In  contrast,  the  clouds  have  less  influence  on  ocean  mixed  layer  during  the 
winter  season  due  to  the  domination  of  shear  production  in  ocean  mixed 
layer  dynamics.  Ihe  effect  of  precipitation  on  ocean  mixed  layer  also  has 
seasonal  variation. 
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Since  both  long  and  short-term  clouds  can  be  shown  to  have  significant 
effect  on  the  dynamic  changes  that  occur  within  the  ocean  mixed  layer,  its 
importance  as  a  central  element  in  any  prognostification  of  the  mixed  layer 
cannot  be  overlooked  or  underestimated. 
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in.  BACKGROUND  DATA 


Ocean  Station  Papa  is  at  the  geographic  position  (50°N,  145°W)  where  a 
Canadian  weather  ship  used  to  be  continuously  stationed.  It  is  located  in  the 
eastern  part  of  the  subarctic  Pacific  region  on  the  southern  edge  of  the  Alaska 
Gyre.  The  dominant  atmospheric  pressure  systems  located  in  the  North 
Pacific  Ocean  are  the  Aleutian  Low  in  winter  and  the  North  Pacific  High  in 
summer.  These  two  systems  control  the  distribution  of  surface  winds  over 
the  Northeast  Pacific  Ocean.  In  winter,  the  Aleutian  Low  directs  air  flow 
northward  into  the  Gulf  of  Alaska,  and  in  summer  the  North  Pacific  High 
directs  winds  southeast  at  Ocean  Station  Papa.  The  direction  of  the  monthly 
winds  does  not  significantly  change  from  season  to  season,  but  the 
magnitudes  of  the  winds  are  about  twice  as  high  in  winter  (13  m/s)  as  in 
summer  (6  m/s).  Because  of  the  large-scale  influence  of  the  wind  systems 
over  the  North  Pacific  Ocean,  the  surface  waters  flow  eastward  across  the 
North  Pacific.  Ocean  Station  Papa  lies  in  the  path  of  this  flow  as  shown  in 
Figure  2. 

According  to  Thompson  (1971),  typical  geostrophic  currents  are  in  an 
east-northeast  direction,  roughly  parallel  to  surface  isotherms  and  isopycnals, 
with  a  speed  around  1  cm/s.  Wind  induced  surface  drift  currents  and  inertial 
currents  may  be  as  large  50  cm/S'  but  since  they  are  typically  either  transient 
or  oscillatory  they  do  not  contribute  appreciably  to  advection  on  a  longer  time 
scale. 
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Figure  2.  Major  Surface  Ocean  Currents  in  the  Sub-Arctic  Region  (adapted 

from  Uda,  1963) 

The  influence  of  horizontal  advection  at  Ocean  Station  Papa  can  bo 
estimated  from  horizontal  currents  and  temperature  gradient.  For  a  typical 
surface  temperature  gradient  in  the  northeast  Pacific  of  about  1°C  (TOO  km)-i 
and  currents  of  1-  10  km  day  ~  the  maximum  temperature  change  to  be 
expected  from  advective  effects  would  be  about  1°C  in  10  days  (Denman  and 
Miyake,  1973).  Tabata  (1965)  attempted  to  estimate  the  temperature  and 
salinity  change  at  Papa  attributable  to  horizontal  advection.  He  concluded 
that  the  average  monthly  temperature  change  at  Papa  due  to  advective  effects 
for  a  five-year  period  was  0.26°C  monlh"^,  with  a  maximum  of  0.78°C 
month"  h  The  average  monthly  salinity  change  was  0.02  ppt  month"  1,  with  a 
maximum  of  0.05  ppt  month"  T  During  the  warming  season,  change  in  sea 
surface  temperature  associated  with  synoptic-scale  weather  patterns  are  of  the 
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order  of  1°C  over  several  days  (Denman  and  Miyate,  1973).  Tabata  (1961,  1965) 
and  Thomson  (1971)  examined  the  vertical  structure  of  the  ocean  at  Papa  in 
some  detail.  At  Papa  the  main  or  permanent  py crocline  occurs  between  100 
and  200  m,  where  the  salinity  increases  from  32.8  ppt  to  33.8  ppt,  the 
temperature  decreases  from  4.5  to  4.0°C,  and  consequently,Ot  increases  from 
26.0  to  26.8.  In  summer,  the  seasonal  thermocline  forms  in  the  upper  75m, 
where  the  temperature  variation  dominates  the  density  variation.  During 
winter  season,  the  salinity  variations  associated  with  the  large  evaporation  at 
the  sea  surface  may  become  significant. 

Tabata  (1961)  indicated  that  there  is  an  annual  average  excess  of 
precipitation  over  evaporation  of  0.5m  in  the  vicinity  of  Ocean  Station  Papa. 
At  Ocean  Station  Papa,  the  monthly  mean  precipitation(P,),  evaporation(E), 
and  E-Pr  are  shown  in  Figure  3.  The  cloud  cover  variation  (about  1/2  during 
winter  and  1  during  summer)  is  shown  in  Figure  4.  These  two  data  sets  are 
the  input  data  used  to  drive  the  NPS  mixed  layer  model. 

There  are  three  reasons  for  using  this  data  set  for  this  study:  (1)  cloudy 
days  occurring  frequently  in  the  vicinity  of  Ocean  Station  Papa,  (2)  salinity 
varies  significantly  in  the  winter  season,  and  (3)  there  are  about  30  years  of 
continuous  meteorological  and  oceanographic  data  in  this  region. 
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Figxire  3.  Monthly  Mean  Values  of  Precipitation,  Calculated  Evaporation,  and 
E  -  Pr  at  Ocean  Station  Papa  (from  Tabata,  1965) 
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Figure  4.  Cloud  Cover  Variation  at  Ocean  Station  Papa  in  1959  (from 

Garwood  and  Adamec,  1982) 
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IV.  BASIC  THEORY  OF  THE  NFS  MIXED  LAYER  MODEL 


A.  OCEAN  MIXED  LAYER  STRUCTURE 

The  upper  layer  of  the  ocean  is  a  region  in  which  the  water  temperature  is 
nearly  isothermal.  In  the  Oceanic  Planetary  Boundary  Layer  (OPBL),  mixing 
due  to  wind  stirring  and  upward  buoyancy  flux  maintains  the  nearly 
homogeneous  temperature  and  salinity  profile.  Depending  on  the  strength  of 
the  winds  and  on  the  direction  of  the  surface  heat  flux,  the  depth  of  the  mixed 
layer  can  be  as  shallow  as  a  few  centimeters  or  as  deep  as  two  hundred  meters 
or  more.  Figure  5  shows  an  idealized  temperature  profile  of  the  upper  ocean. 
An  isothermal  layer  of  temperature  T  exists  in  the  uppermost  h  meters. 
Below  this,  there  is  typically  a  positive  temperature  jump  AT,  over  a  small 
vertical  distance  Ah.  Below  2  =-h-Ah,  the  temperature  decreases  with 
increasing  depth. 


Figure  5.  Idealized  Temperature  Profile  for  Upper  200  m  of  the  Ocean 
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The  characteristics  of  salinity  within  the  mixed  layer  structure  are 
similar  to  those  of  temperature  (Figure  6).  There  is  a  relatively  constant 
salinity  from  z  =  0  to  z  =  -h.  A  corresponding  jump  in  salinity  then  occurs 
across  the  entrainment  zone  from  z  =-h  toz  =  -h-Ah,  followed  by  gradually 
increasing  increments  of  salinity  within  the  depths  below  the  mixed  layer. 
The  sign  of  AS  may  be  positive  or  negative,  depending  upon  the  history  of 
net  precipitation  minus  evaporation,  P,-E. 

z 


Figure  6.  Idealized  Salinity  Profile  for  the  Upper  200  m  of  the  Ocean 

B.  THE  BASIC  EQUATIONS  OF  THE  NPS  MIXED  LAYER  MODEL 

For  most  one-dimensional  ocean  mixed  layer  theories,  the  vertical 
mixing  within  the  turbulent  boundary  layer  and  entrainment  mixing  at  its 
base  occur  in  response  to  local  atmospheric  forces — the  surface  wind  stress 
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and  the  buoyancy  flux  at  the  sea  surface.  The  buoyancy  flux  is  attributable  to 
net  heat  flux  (including  radiation),  evaporation  and  precipitation.  The  shear 
production  of  turbulence  is  attributable  to  surface  wind  stress.  The  mixed 
layer  temperature  Tg,  salinity  Sg,  and  depth  h  are  predicted  by  a  simplified 
form  of  the  NPS  mixed  layer  model  (Garwood,1977)  modified  to  include 
salinity  and  advection: 

dTg  Qo 

=  (2.1) 

-  -U’e  (Ss-S_fc)  +  (E-P,)Ss  (2.2) 


Here  h  is  the  mixed  layer  depth,  Cp  is  specific  heat  for  sea  water  under 
constant  pressure,  p^,  is  the  characteristic  water  density,  E  is  the  surface 
evaporation,  Pr  is  the  precipitation,  and  We  is  the  entrainment  velocity 
parameterized  as; 


^  hg[aAT  -  (5As] 


(2.3) 


where  AT=Tg-T-h  and  AS=Sg-S-h- 

To  solve  equation  (2.3)  for  We,  the  value  of  Iw'l,  T,  h  and  AT  must  be 
known,  together  with  constant  g  and  a  .  From  Garwood  (1977)  the  total 
turbulent  kinetic  energy  (TKE)  is  solved  by  using; 

^(ft{£))  =  G„  +  Gk-Bo-Bj-D.  (2.4) 


where 
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Gn  =  12 


( .  V 


\p 


is  the  wind  stress  production, 

is  the  entrainment  shear  production, 

h  =  «^(Qo/Pft)Cp)-/J^Ss(E-P,) 

is  the  buoyant  damping  or  production  due  to  surface  buoyancy  flux, 

Bh  =  aghATwe 

is  the  buoyant  damping  due  to  entrainment,  and 

_  3 
D  =  2(e)2 

is  dissipation. 

The  time  rate  of  change  in  mixed  layer  depth  is  computed  by 

dh 
dt  ~ 


(2.5) 


(2.6) 


(2.7) 


(2.8) 


(2.9) 


(2.10) 


where  W-h  is  the  upwelling  velocity  specified  at  the  bottom  of  the  mixed 
layer. 


C  CLOUD  EFFECT  ON  HEAT  FLUX  AT  THE  OCEAN  SURFACE 

There  are  two  important  effects  of  clouds  on  the  buoyancy  flux  (Bq)  at  the 
ocean  surface;  (1)  decrease  in  downward  buoyancy  flux  by  reducing  the 
incoming  solar  radiation  or  (2)  decrease  in  upward  buoyancy  fiux  by  emitting 
the  long  wave  radiation  back  to  the  ocean  surface.  Additionally,  the  salinity 
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flux  caused  by  precipitation  also  can  increase  downward  buoyancy  flux 
indirectly. 

The  surface  heat  flux  (upward  positive),  Qo,  is  computed  by 
Qo  =  Qb  ~  Qs  +  ^P<wE  +  Qh 

where  Qs  is  the  incoming  solar  radiation  absorbed  by  the  ocean  surface,  Qt,  is 
the  net  back  radiation,  L  is  the  latent  heat  of  vaporization  of  water,  Qh  is  the 
sensible  heat  flux  to  the  air. 

The  surface  evaporation,  E,  and  sensible  heat  flux  from  the  ocean  surface, 
Qh,  are  estimated  using  the  bulk  aerodynamic  formulae: 

Q/.=P«c|,“’CdU,oCh(T,-T„,)  (2,13) 

where  Cd  is  the  drag  coefficient,  taken  as  0.001  in  this  study,  U]o  is  the  wind 
speed  at  10  m  height,  i}s(T)  is  the  saturated  mixing  ratio,  is  the  air 

{(j) 

temperature  at  the  ocean  surface,  Cp  is  specific  heat  of  the  atmosphere,  Cj],C£ 

are  heat  and  moisture  transfer  coefficients.  In  this  study  we  assume  Cf  =  Cy/. 

Clouds  reduce  the  solar  radiation  upon  the  ocean  surface  by  scattering  and 
absorption  which  is  computed  by 

Qs  =  (l-fla^)(l -1.289  xlO-'^^C'^^jQ,-  (2.14) 

where  C  is  the  proportion  in  eighths  of  a  sky  covered  by  cloud,  and  Qi^^  is  the 
clear  sky  radiation  given  by  (Seckel  and  Beaudry,  1973) 

Qc  A)  ^1  COS0  +  sinp  +  /42COs20  +  Bt  sin2(;f).  (215) 
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The  constant  a  and  b  are  adopted  from  Tabata  (1964)  and  the  cubic  cloud 
cover  correction  from  Laevastu  (1960).  Noon  altitude  of  the  sun  is  a  .  The 
coefficients  (Aq,  Ai,  A^,  B-[,  B2)  were  calculated  by  harmonic  representation  of 


the  values  listed  in  the  Smithsonian  Meteorological  Table,  and 

(  2jt 


(p  = 


365 


(1-21) 


(2.16) 


where  t  is  the  Julian  day  of  the  year. 

The  ocean  surface  emits  longwave  radiation  to  the  atmosphere.  The 
clouds,  as  well  as  dry  air,  partially  absorb  the  radiation  and  reemit  longwave 
radiation  back  to  the  ocean  surface.  Thus  the  net  back  radiation  Q[,,  is 
corrected  for  the  downward  radiation  by  the  clouds.  The  net  back  radiation  is 
estimated  from  the  empirical  formula  (Husby  and  Seckel,  1978) 

Qb  =  1. 14  X 1 0'^  ( 273. 1 6  +  Ts  )'*  (.  39-.  )(l  -  9. 375  x  1 0  -3  C‘ ).  (2.17) 

Here,  Ts  is  sea  surface  temperature,  Ea  is  the  vapor  pressure  of  air  based  on 
the  dew-point  temperature. 

D.  SALINITY  EFFECT  ON  MIXED  LAYER  DYNAMICS 
1.  The  Shallowing  Regime 

When  the  ocean  surface  is  strongly  heated  or  there  is  heavy 
precipitation,  the  mixed  layer  depth  is  taken  as: 

h  =  ^ - .  =  L  (2.18) 

Q  agQo/(p^Cp 
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where  L  is  the  Obokhov  length  scale,  is  the  water  surface  friction  velocity, 

2  2 

which  is  computed  by  u  *  =  (pa/pzi?)Cp  U  If  the  precipitation  and 


evaporation  are  included,  equation  (2.18)  becomes 

Cl 


h  = 


C2asQo/(paCp)-^/3s(P,-E)s' 


(2-19) 


The  heat  equation  for  the  shallowing  regime  is: 

ST  Qo 


(2.20) 


As  a  result,  precipitation  minus  evaporation  influences  mixed  layer 
temperature  indirectly.  For  downward  buoyancy  flux,  a  positive  (Pr-E)  will 
decrease  h,  and  concentrate  the  heat  into  a  shallower  layer.  For  the  case 
where  evaporation  exceeds  precipitation,  a  decrease  of  mixed  layer 
temperature  would  be  expected. 

2.  The  Entrainment  Regime 

For  an  entraining  mixed  layer,  the  buoyancy  discontinuity  at  the 
bottom  of  the  mixed  layer  is: 

Ab  =  agAT  -  PgAS 
the  time  rate  of  change  of  mixed  layer  salinity  is 

OSj  _  Sg(E~P^)  w^ASg 

h 

Equation  (2.22)  demonstrates  that  the  time  rate  of  change  of  salinity  depends 
on  surface  fresh  water  influx  and  entrainment  processes. 


(2.21) 

(2.22) 
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E  METHOD  OF  SOLUTION 

Prediction  of  the  rate  of  deepening  (or  retreat)  of  the  NPS  mixed  layer 
model  depends  upon  an  understanding  of  the  dynamics  of  the  entrainment 
process.  The  turbulence  of  the  overlying  mixed  layer  provides  the  energy 
needed  to  destabilize  and  erode  the  underlying  stable  water  mass  (Garwood, 
1977).  Therefore  the  Turbulent  Kinetic  Energy  (TKE)  budget  is  the  basis  for 
the  entrainment.  This  system  is  closed  using  a  mean-turbuient-field 
modeling  of  the  vertically  integrated  equations  for  the  individual  TKE 
components,  plus  the  inclusion  of  the  bulk  buoyancy  and  momentum 
equations. 

Separate  vertical  and  horizontal  equations  for  TKE  are  used  to  better 
model  the  mixing  process.  Buoyancy  produced  energy  is  somewhat  more 
efficient  than  shear  production  as  a  source  of  energy  for  vertical  mixing 
because  of  its  direct  effect  on  the  vertical  component  of  the  turbulent  velocity. 
The  buoyancy  equauon  is  generated  fr^m  the  heat  and  salt  equations  together 
with  the  equation  of  state: 

p=PQ^l-a(e-0o)  +  p{s-So)\  (2-23) 

and  buoyancy  is  given  by: 

b  =  8{Po-p)/Po  (2-24) 

where  0  is  temperature,  s  is  salinity,  p  is  density,  g  is  gravity,  and  the 
constants  a  and  (3  are  the  expansion  coefficients  for  heat  and  salt.  The  tilae 
represents  the  ‘:otal  instantaneous  value  and  the  subscript  zero  denotes  an 
arbitrary,  but  represert.ativo  constant  value.  Using  b  instead  of  a  0  for  the 
definition  of  buoyancy  allows  this  model  to  be  generalized  to  include  salinity. 
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Temperature  and  salinity  profiles  are  required  as  input  for  model 
initialization.  These  are  used  to  compute  the  mixed  layer  depth,  h.  The 
numerical  solution  for  the  NFS  mdxed  layer  model  defines  the  minimum 
mixed  layer  depth  to  be  1  m.  Other  ocean  environment  parameters  to  be 
specified  include  the  fraction  of  short  wave  radiation  absorbed  in  the  top  one 
meter  of  the  ocean,  the  radiation  extinction  coefficient  for  absorption  of 
radiation  with  depth,  and  the  critical  Richardson  number  for  dynamic 
stability  to  be  maintained  at  the  bottom  of  the  mixed  layer. 

The  parameters  required  for  surface  boundary  condition  computation 
include  wind  speed  and  direction,  cloud  cover,  sea  surface  temperature,  air 
temperature,  dew-point  temperature,  incident  solar  radiation,  and  the  rates  of 
evaporation  and  precipitation.  It  should  be  noted  that  not  all  of  the  incoming 
short  wave  radiation  penetrates  the  ocean  mixed  layer.  Approximately  half 
(for  open  ocean)  is  absorbed  within  the  first  meter.  The  amount  absorbed 
varies  from  region  to  region  depending  on  the  amount  of  absorbing 
particulates  such  as  phytoplankton  and  yellow  substance.  More  radiation  will 
be  absorbed  in  coastal  regions  than  in  the  open  ocean  region.  This  portion  of 
absorbed  short  wave  radiation  is  therefore  considered  to  be  part  of  the  upward 
heat  flux. 

In  a  nondimensional  context,  there  will  be  two  degrees  of  freedom  for  the 

-  3 

dimensionless  entrainment  flux,  =  -  bw  (-h)h/u».  The  independent 

3  - 

parameters  are  H*  =  h/L\  and  Z*  =  h/Lj  where  L]  =  -u*/  bzv  (0)  (the  Obokhov 

length  scale)  and  L2  =  «*//(/  is  coriolis  parameter).  New  ocean  mixed  layer 
depth,  temperature  and  salinity  are  predicted  at  one  hour  interval.'^.  The  steps 
in  the  prediction  computation  are  shown  in  Figure  7. 
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V.  PROCEDURES,  RESULTS  AND  DISCUSSION 


A.  PROCEDURE 

The  aim  of  this  study  is  to  examine  the  effects  of  cloud  cover  and 
precipitation  on  the  ocean  mixed  layer.  Data  collected  in  1959  was  used  in  this 
study.  Three  experiments  with  different  time  scales  are  designed  for  this 
study:  namely,  short-term  experiment  (3-day  period),  medium-term 
experiment  (30-day  period)  and  long-term  experiment  (365-day  period).  Each 
experiment  is  composed  of  three  different  cases  (the  standard  case.  Case  1  and 
Case  2). 

In  the  short-term  experiment,  the  standard  case  is  as  follows:  the  cloud 
cover  is  daily  mean  values  of  the  three-hourly  observational  data;  the 
precipitation  minus  evaporation  rate  (E-P,)  is  assumed  to  be  0.25  cm/day.  The 
case  with  20%  reduction  of  cloud  cover  to  the  standard  case  with  the  same 
rate  of  precipitation  minus  evaporation  is  taken  as  Case  1.  The  case  with  the 
same  amount  of  cloud  cover  as  the  standard  case  but  with  E-P^  =  0  is  taken  as 
Case  2  (Table  1).  In  this  experiment  we  focus  on  the  cloud  and  precipitation 
effect  on  the  diurnal  variation  of  the  mixed  layer  depth  (MED)  and 
temperature  (Ts). 
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TABLE  1.  SIMULATION  CASES  FOR  THE  SHORT-TERM  EXPERIMENTS 


Month 

Days 

Standard  Case 

Case  1 

Case  2 

Cloud 

Cloud 

Qoud 

Cover 

Cover 

Cover 

(1/8) 

(1/8) 

(1/8) 

25 

ra 

0.25 

6 

0.25 

ra 

0 

Jan. 

26 

5.3 

0.25 

4.2 

0.25 

5.3 

0 

27 

JBBk 

0.25 

5.8 

0.25 

ra 

0 

156 

8 

0.25 

ra 

0.25 

8 

0 

Jun. 

157 

8 

0.25 

19 

0.25 

8 

0 

158 

8 

0.25 

0.25 

8 

0 

As  for  the  short-term  experiment,  three  cases  are  designed  for  the 
medium-term  experiment,  where  daily  mean  value  of  cloud  cover  and 
precipitation  are  replaced  by  monthly  mean  value  in  the  standard  case  (Table 
2).  The  monthly  mean  value  of  cloud  cover  comes  from  three-hourly 
observations  while  the  montly  mean  value  of  precipitation  are  specified  by 
Figure  3. 


TABLE  2.  SIMULATION  CASES  FOR  THE  MEDIUM-TERM  EXPERIMENTS 


Month 

Days 

Standard  Case 

Case  1 

Case  2 

Qoud 

■ra 

Qoud 

Pr 

Ev 

Cloud 

Prtv 

Cover 

^^SJra 

Cover 

(cm/ 

(cm/ 

Cover 

(cm/ 

(1/8) 

(1/8) 

month) 

month) 

(1/8) 

month) 

Jan. 

1-31 

6.3 

12.5 

2.5 

5.0 

12.5 

2.5 

6.3 

0 

Jun. 

151- 

ra 

2.8 

1.25 

5.9 

2.8 

1.25 

5.9 

0 

181 

ra 

Two  months  that  could  be  characterized  as  typical,  nonetheless  produced 
interesting  results.  The  first  is  June,  where  a  strong  net  heat  gain  (Bo  >  0)  as 
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well  as  a  weak  wind  speed  (Figure  20)  was  noted.  Also  in  January  a  significant 
net  heat  loss  (Bq  <  0)  and  stronger  wind  speed  were  observed  (Figure  16). 

As  for  the  long-term  (annual)  experiment,  the  time  series  of  changes  in 
MLD  and  Tg  was  composed  with  cloud  cover  and  precipitation  variations. 
Fiere  a  complex  event  with  monthly  mean  value  of  cloud  cover  and 
precipitation  (Table  3)  will  be  compared  with  a  simple  event  using  mean 
yearly  value  of  cloud  cover  (n  =  6.7),  precipitation  (P,  =  62.2  cm/year),  and 
evaporation(Ev  =  48.7  cm/year).  Both  complex  and  simple  events  will  be 
compared  with  actual  bathythermogram(BT)  observations.  In  addition  to 
cloud  cover  and  precipitation,  the  surface  boundary  condition  including  wind 
direction,  wind  speed,  sea  surface  temperature,  air  temperature  and  dew¬ 
point  temperature  are  taken  from  observational  data,  which  were  observed 
every  three  hours.  The  heat  and  momentum  fluxes  are  computed  from  these 
actual  observations.  The  initial  temperature  profile  is  selected  from  one  of 
665  BT  measurements.  This  temperature  profile  plus  observed  precipitation 
and  calculated  evaporation  were  used  as  initial  condition  of  the  model.  In 
these  experiments,  the  sta’idard  case  (subscript  0)  will  be  compared  with  Case 
1  (subscript  1)  and  Case  2  (subscript  2).  When  subtracting  the  results  of  Case  1 

and  Case  2  by  standard  case,  the  differential  values  provide: 

1.  Differential  MLD  (Ho -Hi)  and  differential  Tg  (Tgo-Tgi)  in  variable 
cloud  simulation. 

2.  Differential  MLD  (Hq  -  H2),  differential  Tg  (Tgo  -  Tg2)  and  differential 
salinity  (So  -  S2)  in  variable  precipitation  simulation. 

The  working  hypothesis  here  is  that  cloud  cover  change  will  make  sea 
surface  temperature  change,  which  will  lead  to  mixed  layer  and  boundary 
condition  changes.  Thus,  the  MLD  and  Tg  prediction  will  become  more 
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complicated.  In  order  to  simplify  the  problem,  we  assume  all  boundary 
variables  were  constant  except  cloud  cover  and  precipitation  in  these 
experiments. 


TABLE  3.  SIMULATION  CASES  FOR  THE  LONG-TERM  EXPERIMENT 

(COMPLEX  EVENT) 


Month 

Days 

Standard  Case 

Case  1 

Case  2 

Ooud 

Cloud 

Pr 

Ev 

Qoud 

PrEv 

Cover 

9R^9 

Cover 

(cm/ 

(cm/ 

Cover 

(cm/ 

(1/8) 

(1/8) 

month) 

month) 

(1/8) 

month) 

Jan. 

1-31 

6.3 

12.5 

2.5 

5.0 

12.5 

2.5 

6.3 

0 

Feb. 

32-59 

5.9 

2.5 

4.8 

IB 

2.5 

4.8 

5.9 

0 

Mar. 

60-90 

6.8 

3.8 

wM 

m 

3.8 

6.8 

0 

^9 

91-120 

6.5 

2.5 

1.3 

5.2 

1.3 

6.5 

0 

121-151 

^9 

2.5 

1.9 

5.8 

2.4 

1.9 

^9 

0 

Jun. 

152-181 

2.8 

1.25 

5.9 

2.8 

1.25 

7.4 

0 

Jul. 

182-212 

m 

5.0 

1.2 

6.1 

5.0 

1.2 

7.6 

0 

^9 

213-243 

m 

2.5 

5.9 

4.4 

2.5 

0 

Sep. 

244-273 

6.6 

2.5 

3.0 

5.3 

2.5 

3.0 

6.6 

0 

Oct. 

274-304 

6.3 

8.8 

5.0 

5.0 

8.8 

5.0 

6.3 

0 

Nov. 

305-334 

6.2 

5.6 

6.9 

6.9 

5.6 

6.9 

6.2 

0 

Dec. 

335-365 

5.8 

5.6 

15.0 

4.6 

5.6 

15.0 

5.8 

0 

B.  RESULTS  AND  DISCUSSION 
1.  Short-Term  Experiments 

Two  events  were  extracted  from  the  months  of  January  and  June. 
These  events  were  chosen  at  random,  and  as  such,  do  not  represent 
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exceptional  or  special  cases.  Diurnal  variation  of  MLD,  Tj  and  mixed  layer 
salinity  were  examined  in  these  periods. 

a.  Event  #1  (January,  days  15-27) 

The  wind  speed  was  low  at  the  beginning  of  the  event  and  later 
reached  its  maximum  speed  in  the  middle  of  the  second  day  (Figure  8a).  The 
model  simulation  began  an  hour  before  sunset.  The  initial  MLD  was  not 
deep  (45  m)  when  compared  to  the  typical  MLD  in  January.  The  net  surface 
heat  flux  indicated  a  net  heat  loss  during  this  period  (Figure  8b).  Figure  9 
shows  the  variation  of  differential  MLD  and  Tg  in  the  variable  cloud 
simulation.  The  differential  MLD  was  negative  during  each  night  and 
positive  during  each  daylight  period.  The  maximum  negative  differential 
MLD  occurs  after  midnight  of  the  first  day  due  to  a  lull  in  the  wind  velocity. 
A  near  zero  differential  MLD  appeared  the  second  day  in  association  with  the 
onset  of  stronger  winds.  The  key  to  examining  this  phenomenon  is  the 
disproportionate  darkness  of  the  long  January  nights  (~16  hours).  The 
greenhouse  effect  was  particularly  active  during  the  lengthy  hours  of 
darkness,  resulting  in  pronounced  long-wave  radiation  return  back  to  the 
ocean  surface  in  the  standard  case,  thereby  producing  a  shallower  MLD. 

After  the  sun  rise,  the  albedo  effect  dominated  variations  of  the 
MLD  when  the  surface  buoyancy  flux  was  downward.  Then,  due  to  a  net 
increase  in  solar  radiation  reaching  the  ocean  surface  in  Case  1,  the 
interaction  of  these  dynamics  leads  to  a  shallower  MLD  than  would  be  found 
in  the  standard  case. 

The  significant  greenhouse  effect  throughout  the  hours  of 
darkness  will  result  in  a  warming  of  the  ocean  surface,  a  form  of  positive 
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differential  Tj.  Negative  differential  Ts,  found  during  daylight  hours,  were 
associated  with  the  deeper  MLD  of  the  standard  case.  The  maximum 
differential  of  T*  occurred  toward  the  end  of  daylight  on  the  first  day  when 
the  MLD  was  shallow. 

Unlike  the  variations  in  differential  MLD  and  Tg  found  in  the 
variable  cloud  simulation,  the  differential  MLD  and  Tg  in  the  variable 
precipitation  simulation  produced  consistently  negative  values  throughout 
the  period  (Figure  10).  The  influx  of  fresh  water  into  the  mixed  layer  will 
result  in  a  shallower  MLD,  ultimately  producing  a  lower  mixed  layer 
temperature  associated  with  a  net  heat  loss  in  this  period.  Again,  the 
maximum  negative  differential  MLD  and  Tj  noted  during  the  variable 
precipitation  simulation  are  associated  with  a  lull  in  the  wind  occurring 
during  the  hours  of  darkness  on  the  first  day.  The  differential  of  mixed  layer 
salinity  increased  in  a  stepwise  manner  throughout  the  period  (Figure  lla).  It 
should  be  noted  that  variations  in  differential  salinity  are  dependent  upon 
variable  wind  velocities,  i.e.  in  the  case  of  E-P,  <  0,  an  increase  of  wind 
velocity  will  entrain  saline  water  from  the  lower  layers  into  the  mixed  layer, 
thereby  increasing  the  mixed  layer  salinity  and  decreasing  the  differential 
salinity. 

Comparing  the  mean  value  of  MLD  and  T^  in  this  event  (see 
Table  4),  two  features  arise  worthy  of  note.  For  the  first,  when  we  compare 
the  difference  in  mean  MLD  between  the  standard  case  and  Case  1  (-0.3  m), 
the  difference  is  strikingly  less  than  what  appears  in  a  comparison  of  the 
standard  case  and  Case  2,  when  a  differential  of  about  4  m  is  found.  The 
implication  is  that  precipitation  has  had  a  significant  and  pronounced  effect 
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upon  the  MLD  during  this  event.  Secondly,  we  may  note  that  due  to  a  deeper 
MLD  and  weaker  solar  radiation,  the  mean  value  of  Tj  has  nearly  no  change 
throughout  the  period. 


TABLE  4.  THREE-DAY  MEAN  VALUES  OF  MLD,  TEMPERATURE,  NET 
SURFACE  HEAT  FLUX  AND  MIXED  LAYER  SALINITY 


Event 

Days 

Case 

IQSI 

25 

Standard 

45.5 

5.4 

33.8 

31.9 

January 

1 

1 

45.9 

5.4 

32.9 

31.9 

27 

2 

49.6 

5.4 

33.8 

32.0 

156 

Standard 

36.1 

7.6 

-73.4 

31.9 

June 

/ 

1 

32.4 

7.5 

-160.9 

31.9 

158 

2 

37.1 

7.6 

-73.4 

32.0 

b.  Event  #2  (June,  days  156-158) 

In  this  event  the  MLD  is  approximately  40  m,  and  the  downward 
buoyancy  flux  was  stronger  than  event  #1.  The  wind  was  more  pronounced 
at  the  beginning  of  the  event  but  subsequently  decreased  from  the  middle  of 
the  first  day  (Figure  12).  Daylight  and  the  associated  period  of  insolation 
comprised  13  hours.  The  model  simulation  began  in  the  afternoon,  four 
hours  before  sunset. 

The  differential  MLD  and  Tj  in  the  variable  cloud  simulation  are 
displayed  in  Figure  13.  At  the  beginning  of  the  event,  a  stronger  wind  and 
weak  upward  buoyancy  flux  lead  to  a  near-zero  differential  MLD  and  Tj 
during  the  hours  of  darkness  of  the  first  day.  Then  as  conditions  change,  the 
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positive  differential  MLD  and  the  negative  differential  Tg  increased  during 
the  day  and  diminish  in  the  darkness.  This  pattern  indicates  that  the  albedo 
effect  has  significant  influence  on  MLD  and  Tg  during  any  day  having  strong 
solar  radiation  and  a  lengthened  duration  of  insolation.  In  the  standard  case, 
the  clouds  scatter  and  absorb  a  significant  amount  of  incoming  solar 
radiation,  decreasing  the  downward  buoyancy  flux  at  the  ocean  surface  during 
the  day.  This  decrease  in  the  downward  buoyancy  flux  will,  in  turn,  lead  to 
an  increase  in  entrainment  velocity  thereby  deepening  the  mixed  layer  and 
cooling  the  ocean  surface.  The  greenhouse  effect  had  little  significant  impact 
on  the  mixed  layer  during  the  period  of  this  event  due  to  relative  decrease  in 
the  hours  of  darkness  and  weaker  long  wave  radiation. 

The  differential  MLD  produced  negative  values  in  the  variable 
precipitation  simulati^'n  higure  14a).  The  precipitated  fresh  water  has  the 
effect  of  increasir^  ratification  of  the  upper  ocean  in  the  standard  case.  As 
the  stratification  becomes  more  pronounced,  so  does  the  density  jump  at  the 
base  of  the  mixed  layer,  resulting  in  a  decrease  in  the  rate  of  entrainment. 
With  the  decrease  in  the  rate  of  entrainment,  the  MLD  then  retreated  relative 
to  the  case  of  E-P,  -  0.  A  positive  differential  Ts  is  associated  with  a  shallower 
MLD  in  the  standard  case  (Figure  14b).  The  variation  in  differential  salinity 
grew  in  stair-step  fashion  throughout  the  period  of  the  event.  (Figure  lib). 
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Figure  8.  Variation  of  (a)  Wind  Stress  and  (b)  Net  Surface  Heat  Flux  (January, 

days  25-27) 
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Figure  9.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  Ts  in  the 
Variable  Cloud  Simulation  (January,  days  25-27) 
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Figure  10.  Varialion  of  (a)  Differential  MLD  and  (b)  Differential  in  the 
Variable  Precipitation  Simulation  (January,  days  25-27) 
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Figure  11.  Variation  of  Differential  Salinity  in  (a)  Event  #1  and  (b)  Event  #2 

for  the  Short-term  Experiment 
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Figure  12.  Variation  of  (a)  Wind  Stress  and  (b)  Net  Surface  Heat  Flux  (June 

days  156-158) 
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Figure  I3.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  Tj  in  the 
Variable  Cloud  Simulation  (June,  days  156-158) 
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Figure  14.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  in  the 
Variable  Precipitation  Simulation  (June,  days  156-158) 
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Examining  the  mean  value  of  the  net  surface  heat  flux  (see  Table 
4),  we  find  a  significant  difference  between  the  standard  case  and  Case  1 
during  the  event.  The  mean  surface  heat  flux  shows  a  large  increase,  87.5 
w/m2  in  this  event,  compared  to  0.9  w/m^  in  Event  #1.  This  indicates  that 
the  clouds  have  significant  impact  on  surface  buoyancy  flux  and  directly 
influence  upper  ocean  dynamics  during  this  event.  Another  striking  feature 
that  arises  from  the  figures  in  Table  4  is  the  difference  of  mean  MLD  between 
the  standard  case  and  Case  2  (0.9m)  that  is  distinctly  less  than  those  of  Event 
#1  (4.008m).  These  figures  illustrate  that  the  dynamic  effect  of  precipitation 
on  MLD  was  greater  for  the  circumstances  of  Event  #1. 

2.  Medium-term  Experiments 
a.  January  (days  1-31) 

In  general,  the  MLD  is  deeper  than  the  annual  average  in  this 
period,  subject  to  interaction  with  strong  winds  and  strong  upward  buoyancy 
flux  (Figures  15a  and  16).  The  model  simulation  began  at  an  hour  before 
dark. 

The  variation  of  differential  MLD  in  the  variable  cloud 
simulation  produced  a  positive  differential  MLD  during  the  day  and  a 
negative  differential  MLD  during  the  night  (Figure  17a).  This  result  is  similar 
to  the  results  found  in  the  short-term  experiment  of  Event  #1,  where  the 
greenhouse  effect  is  the  dominant  factor  on  the  MLD  during  the  night  while 
the  albedo  effect  became  the  significant  factor  during  the  day,  associated  with  a 
downward  buoyancy  flux.  A  larger  differential  MLD  than  what  might 
otherwise  be  expected  is  a  reflection  of  a  lull  wind  which  can  serve  to  either 
deepen  or  shallow  the  MLD. 
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With  the  exception  of  a  few  negative  differential  Tj's  occurring  in 
daytime,  and  which  are  essentially  anomalies,  the  differential  was 
generally  positive  throughout  the  period  (Figure  17b).  Varying  consistently 
with  the  diurnal  cycle,  decreasing  differential  Ts  correlated  with  daylight. 
Differential  Ts  increased  consistently  with  night.  The  differential  Ts  in  this 
circumstance  is  a  clear  reflection  of  the  predominant  greenhouse  effect  in 
January.  The  stronger  upward  buoyancy  flux  associated  with  the  mid-latitude 
winter  regime  will  result  in  a  colder  ocean  surface  in  Case  1  as  more  long 
wave  radiation  escapes  from  the  ocean  surface.  In  the  case  of  shallower  MLD 
in  daylight  hours  due  to  weaker  wind  velocity,  the  ocean  surface  in  Case  1 
may  become  warmer  and  form  a  negative  differential  Tg. 

The  variation  of  differential  MLD  and  T^  in  the  variable 
precipitation  simulation  are  displayed  in  Figure  18.  As  a  result  of 
precipitation  damping  turbulent  kinetic  energy  and  consequently  inducing 
the  MLD  to  become  shallower  in  the  standard  case,  a  negative  differential 
MLD  appears  throughout  the  period.  The  shallower  MLD  in  the  standard 
case  then  lost  more  upward  long  wave  radiation,  and  caused  ocean  surface  to 
be  colder  than  is  the  case  with  E-P,  =  0.  Quite  predictably,  a  negative 
differential  Tj  also  appeared  throughout  the  period.  There  was  one  exception 
to  this  otherwise  rigid  pattern,  however.  During  the  daylight  hours  of  the  day 
16,  a  single  instance  of  positive  differential  Tj  occurred  due  to  the  shallowest 
MLD  that  developed  during  the  period,  which  increases  the  mixed  layer 
temperature  in  the  standard  case  A  striking  variation  in  differential  salinity 
also  occurred  on  day  16  while  the  MLD  was  at  its  shallowest  (Figure  19a). 
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b.  June  (Days  152-181) 


Except  for  a  storm  that  occurred  on  the  second  day  of  the  month, 
the  winds  were  moderate  until  the  middle  of  the  month  and  minimal 
afterwards  (Figure  20a).  The  value  of  net  surface  heat  flux  indicate  a  strong 
downward  buoyancy  flux  at  the  ocean  surface  during  the  period  (Figure  20b). 
The  MLD  was  approximately  35m  at  the  beginning  of  the  period  and 
correspondingly  decreased  after  day  13  as  wind  velocities  began  to  subside 
(Figure  15b).  The  model  simulation  began  in  the  afternoon,  four  hours  before 
darkness. 

Figure  21  shows  the  variation  of  differential  MLD  and  T^  in  the 
variable  cloud  simulation.  Prior  to  day  13,  a  few  hours  of  negative 
differential  MLD  occurred  in  the  hours  of  darkness.  But  after  day  13,  the 
differential  MLD  became  consistently  positive  and  established  a  pattern  of 
variation  in  strict  conformance  with  variations  in  wind  magnitude,  higher 
wind  velocities  increasing  the  differential  MLD.  This  indicates  that  the 
albedo  effect  dominates  the  variation  of  MLD  in  June,  associated  with  a 
typical  mid-latitude  summer  regime  (i.e.  strong  solar  radiation  and  longer 
insolation  duration). 

The  negative  differential  T^  also  kept  a  constant  value  before  day 
13  when  the  MLD  was  relatively  deeper.  The  maximum  differential  Tg  (in 
excess  of  4°C)  is  reached  on  the  day  27  when  the  MLD  becomes  quite  shallow 
(<5m).  This  sequence  suggests  that  the  cloud  effect  on  mixed  layer 
temperature  is  particularly  dramatic  when  the  mixed  layer  is  very  shallow. 
After  day  27,  the  MLD  deepens  again  due  to  increasing  wind  velocity,  and  the 
differential  T^  therefore  decreased  to  near  zero. 
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The  variation  of  differential  MLD  and  Tg  in  the  variable 
precipitation  simulation  is  shown  in  Figure  22.  Similar  to  result  showed  in 
January,  the  negative  differential  MLD  appeared  throughout  the  period,  but 
the  magnitude  of  the  differential  was  less  than  what  occurred  in  January  due 
to  a  smaller  quantity  of  precipitation  in  June  (see  Table  2).  Associated  with 
the  negative  differential  MLD  there  is  a  positive  differential  Tj  which  holds  a 
constant  value  close  to  zero  preceding  day  13,  and  then  rapid  increase  on  the 
seventeenth  and  twenty-ninth  days  when  the  MLD  becomes  very  shallow.  A 
similar  pattern  of  occurrence  also  has  been  found  in  the  differential  salinity 
variation  (Figure  19b).  This  is  consistent  with  the  fact  that  a  shallow  MLD 
may  lead  to  an  increase  of  variability  for  both  temperature  and  salinity  within 
the  mixed  layer. 

The  mean  value  of  MLD,  T*,  net  surface  heat  flux  and  mixed  layer 
salinity  for  January  and  June  are  shown  in  Table  5.  These  values  would 
indicate  that  the  effects  of  the  clouds  and  precipitation  will  result  in 
interactive  dynamics  on  the  mixed  layer  that  differ  with  the  seasons.  First, 
the  difference  in  mean  MLD  between  the  standard  case  and  Case  1  in  June 
(2.34m)  is  greater  than  that  of  January  (0.56m),  but  the  difference  in  mean 
MLD  between  the  standard  case  and  Case  2  in  June  (0.26m)  is  significantly  less 
than  that  found  in  January  (7.39m).  The  demonstrated  variation  in  these 
results  suggest  that  the  clouds  have  the  greatest  impact  on  MLD  in  June, 
while  the  effects  of  precipitation  were  obviously  more  significant  in  January. 
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Figure  16.  Variation  of  (a)  Wind  Stress  (b)  Net  Surface  Heat  Flux  in  January 
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Figure  17.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  T*  in  the 
Variable  Cloud  Simulation  in  January 
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Figure  IS.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  T*  in  the 
Variable  Precipitation  Simulation  in  January 
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Figure  19.  Variation  of  Differential  Salinity  in  (a)  January  and  (b)  June 
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Figure  21.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  Tj  in  the 
Variable  Cloud  Simulation  in  June 
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Figure  22.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  Tj  in  the 
Variable  Precipitation  Simulation  in  June 
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TABLE  5.  MONTHLY  MEAN  VALUES  OF  MLD,  TEMPERATURE,  NET 
SURFACE  HEAT  FLUX  AND  MIXED  LAYER  SALINITY 


Event 

Days 

Case 

HBl 

1 

Standard 

49.0 

5.9 

49.2 

31.9 

January 

1 

1 

49.6 

5.9 

51.2 

31.9 

31 

2 

56.4 

5.9 

49.2 

32.0 

152 

Standard 

17.1 

9.2 

-109.4 

31.9 

June 

1 

1 

14.7 

10.5 

-179.3 

31.9 

181 

2 

17.3 

9.2 

-109.4 

32.0 

Second,  the  positive  mean  net  surface  heat  flux  indicates  a  net  heat  loss  at  the 
ocean  surface  in  January,  while  a  negative  mean  net  surface  heat  flux  in  June 
denotes  a  net  heat  gain.  A  pronounced  net  heat  loss  will  cause  the  MLD  of 
Case  1  to  be  deeper  than  the  MLD  of  standard  case  due  to  the  lessening  of  the 
greenhouse  effect.  A  reverse  effect  occurs  as  significant  net  heat  gain  leads  to 
a  shallower  MLD  in  Case  1  than  in  the  standard  case  as  the  albedo  effect 
becomes  enhanced  in  the  standard  case.  Another  meaningful  index  of 
seasonal  tendencies  is  the  very  large  difference  between  the  mean  net  surface 
heat  flux  recorded  in  June  (69.8  w/m^)  and  the  substantially  smaller  figures  of 
January  (2.1  w/m^).  The  implication  of  this  large  difference  is  that  variations 
of  the  clouds  in  June  dominate  the  strength  of  surface  buoyancy  flux  and 
directly  influences  mixed  layer  dynamics. 

3.  Long-Term  Experiments 

A  mixed  layer  simulation  of  365  days  was  conducted  in  order  to 
examine  the  cumulative  and  serial  effects  of  cloud  cover  and  rainfall  events 
over  time.  Monthly  mean  value  of  cloud  coverage  and  precipitation  are 
available  for  what  we  have  termed  the  "complex"  event,  and  the  annual 
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mean  value  of  cloud  coverage  and  precipitation  are  used  for  a  "simple" 
event.  Comparing  these  two  events,  we  find  some  distinctions  in  the  MLD, 
temperature  and  salinity  prediction  between  the  two  events.  The  output  of 
the  MLD  and  temperature  from  the  complex  and  simple  events  are  compared 
to  BT  observations. 

a.  Complex  Event 

During  the  annual  cycle  of  the  sea  surface  temperature,  the  ocean 
surface  warms  from  April  through  September  and  cools  during  the 
remaining  months.  The  time  series  of  yearly  variations  in  net  heat  flux 
(Figure  23a)  corresponds  to  this  annual  cycle.  Also  associated  with  the  annual 
cycle  there  is  a  pronounced  downward  buoyancy  flux  in  the  warming  season, 
followed  by  a  shift  to  an  upward  buoyancy  flux  as  the  ocean  surface  cools.  The 
annual  variation  of  surface  wind  stress  is  shown  in  Figure  24.  Cloud  cover  is 
relatively  diminished  during  the  cooling  season,  while  cloudy  days  frequently 
appet  r  in  the  warm  season  (Figure  4).  These  atmospheric  boundary 
conditions  are  in  concourse  with  a  deepening  of  MLD  and  lower  mixed  layer 
temperature  during  the  cooling  season  and  a  shallower  MLD  and  higher 
mixed  layer  temperature  during  the  warming  season  (see  dotted  line  in 
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Figure  23.  Variation  of  Net  Surface  Heat  Flux  in  (a)  Complex  Event  and  (b) 

Simple  Event  in  the  Standard  Case 
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Figure  24.  Variation  of  Wind  Stress  in  1959 

Figure  30).  The  annual  variation  of  differential  MLD  and  in  the  variable 
cloud  simulation  are  shown  in  Figure  25.  An  examination  of  this  data  will 
further  illustrate  some  of  the  difference  between  the  cooling  and  warming 
season.  Firstly  we  may  note  that  negative  differential  MLD  frequently  appears 
in  the  cooling  season,  while  the  warming  season  invariably  produces  a 
positive  differential  MLD.  Secondly,  the  differential  Ts  decreases  from 
October,  the  month  in  which  the  ocean  surface  begins  to  cool,  and  maintains 
a  constant  value,  near  zero,  throughout  the  winter.  In  contrast,  the 
differential  Ts  begins  to  increase  from  the  spring  transition  and  reaches  its 
maximum  value  (~7.0°C)  in  the  late  summer.  Taken  together,  these  two 
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trends  suggest  that  the  albedo  effect  dominates  the  mixea  layer  dynamics 
during  the  warming  season  while  the  greenhouse  effect  plays  an  important 
role  on  mixed  layer  variation  during  the  cooling  season. 

As  summer  approaches,  the  MLD  retreats  in  association  with  an 
increase  in  downward  buoyancy  flux  and  reduced  winds.  For  the  shallower 
MLD  in  Case  1,  the  solar  energy  is  concentrated  in  the  thinner  layer  and  thus 
produces  a  greater  negative  differential  T*. 

During  the  spring  transition  (i.e.  period  between  day  100  and  day 
120)  the  differential  MLD  is  nearly  zero.  In  contrast,  the  differential  Ts 
abruptly  increases.  A  similar  circumstance  develops  in  late  spring  between 
day  177  and  day  180.  The  retreat  of  the  MLD  to  a  thin  layer  at  the  ocean 
surface  was  true  for  both  the  standard  case  and  Case  1  during  the  spring 
transition.  The  unusually  small  MLD  will  uniformly  result  in  a  lesser 
differential  MLD  value  while  the  differential  Ts  will  increase.  After  late 
summer  the  MLD  began  to  deepen  and  a  larger  differential  appeared,  while 
the  differential  Tj  decreased  with  the  deeper  MLD. 
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Figure  25.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  T*  in  the 
Variable  Cloud  Simulation  (Complex  Event) 
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Figure  26.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  T*  in  the 
Variable  Precipitation  Simulation  (Complex  Event) 
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Figure  27.  Salinity  Variation  in  Standard  Case:  The  dotted  line  indicates  the 
complex  event;  the  solid  line  indicates  the  simple  event 

The  annual  variation  of  differential  MLD  and  Ts  for  the  variable 
precipitation  experiments  is  displayed  in  Figure  26.  A  positive  differential 
MLD  appeared  during  the  cooling  season  (February,  November,  and 
December)  when  the  rate  of  evaporation  was  greater  than  precipitation.  The 
largest  negative  differential  MLD  occurs  during  the  months  of  January,  March 
and  October  when  the  precipitation  is  greatest  (see  Table  3).  The  decline  of 
precipitation  during  the  warming  season  results  in  a  decrease  of  differential 
MLD,  while  precipitation  becomes  important  in  its  effect  on  the  mixed  layer 
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from  late  autumn,  increasing  differential  MLD  and  giving  rise  to  an  annual 
cycle  of  variation  in  the  differential  MLD. 

During  the  winter  months,  a  condition  of  negative  differential  Tj 
prevails  as  a  result  of  precipitation  exceeding  evaporation.  This  induces  the 
MLD  in  the  standard  case  to  be  shallower  than  that  found  in  the  case  with  E- 
Pr  =  0,  thereby  enhancing  the  rate  of  cooling  associated  with  winter  time  net 
surface  heat  loss.  Note  that  the  magnitude  of  differential  Tg  in  the  variable 
precipitation  experiment  is  larger  than  that  for  the  variable  cloud  experiment 
during  the  spring  transition.  The  explanation  is  related  to  the  fact  that  solar 
radiation  is  still  relatively  weak  in  the  early  spring. 

The  variation  of  salinity  in  Case  2  largely  depend  upon  variation 
in  wind  velocity  (see  dotted  line  in  Figure  27),  with  large  wind  increasing  the 
mixed  layer  salinity  in  the  case  with  E-P,  <  0.  This  is  exemplified  by  an  abrupt 
increase  in  mixed  layer  salinity  that  occurs  at  the  beginning  of  March  (day  60) 
when  wind  velocity  is  at  its  strongest  for  the  period  (Figure  24).  Similar  to 
variations  in  T*,  mixed  layer  salinity  decreases  at  a  more  rapid  rate  during  the 
spring  transition  and  summer  months  when  wind  velocities  are  low  and  the 
MLD  is  shallow.  When  late  autumn  approaches,  the  mixed  layer  salinity 
tends  to  increase  due  to  erosion  of  seasonal  halocline. 

b.  Simple  Event  (n  =  6.7,  Pr  =  62.7  cm/year,  Ey  =  4S.7  cm/year) 

Comparing  the  net  heat  flux  between  the  complex  event  and  this 
simple  event  (Figure  23),  we  note  that  during  the  summer  the  downward 
buoyancy  flux  in  the  simple  event  is  more  pronounced  than  in  the  complex 
event.  In  the  variable  cloud  simulation  the  results  are  similar  to  those  of  the 
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complex  event,  with  the  exception  of  smaller  magnitudes  of  differential  MLD 
and  Ts  (Figure  28). 

The  variable  precipitation  simulation,  however,  causes  a 
sequence  of  features  quite  unlike  those  found  in  the  complex  event  (Figure 

29) .  First,  the  differential  MLD  has  a  negative  value  throughout  the  entire 
year,  whereas  a  positive  differential  MLD  occurs  in  the  complex  event  during 
the  months  with  an  E-P^  >  0  configuration.  Second,  during  the  late  autumn, 
the  decreasing  rate  of  differential  Tg  is  smaller  in  the  simple  event.  Third,  in 
the  simple  event,  the  mixed  layer  salinity  decreases  only  slightly  during  the 
winter  and  then  abruptly  increases  on  day  60  with  the  arrival  of  high  wind 
speed.  Afterward,  mixed  layer  salinity  falls  into  a  pattern  of  rigid 
correspondence  with  the  variations  of  wind  velocity,  reaching  a  constant 
value  in  late  autumn  (Figure  27).  From  Figure  27,  we  find  that  the  model 
prediction  for  the  complex  event  yields  at  a  much  closer  approximation  of  the 
actual  salinity  profile  (i.e.  large  evaporation  in  cooling  season)  than  does  the 
prediction  for  the  simple  event  . 

Drawing  from  BT  observations,  the  actual  values  of  MLD  and  Tg 
were  compared  to  those  values  predicted  by  the  NFS  mixed  layer  model  for 
both  the  complex  and  simple  events  using  the  standard  case  category  (Figure 

30) .  A  distinction  exists,  in  that  the  BT  data  shows  a  deeper  MLD  profile 
during  the  winter  months  for  both  simple  and  complex  events  and  close 
agreement  with  the  model  output  for  both  events  during  the  summer 
months. 

In  the  general  scheme,  the  values  of  the  model  mixed  layer 
temperature  compare  favorably  with  observed  data.  During  the  winter 
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months,  the  BT  observations  reveal  lower  temperatures  than  the  model 
predictions  in  both  complex  and  simple  events.  And  as  we  progress  through 
the  study  period,  the  actual  temperatures  increase  more  than  the  model 
predictions  for  the  complex  event,  which  is  commensurate  with  the  actual 
MLD  being  shallower  than  the  model  prediction.  During  the  same  time 
period,  circumstances  developed  wherein  a  constant  cloud  cover  and 
precipitation  resulted  in  an  extremely  shallow  MLD,  yielding  the  highest 
temperatures  of  summer.  With  the  approach  of  winter,  the  actual 
temperatures  are  again  lower  than  those  of  the  model  prediction  for  both  the 
complex  and  the  simple  events. 

A  second  comparison  between  the  actual  recorded  data  and  the 
model  prediction  can  be  made  by  examining  the  representations  of  mean 
value  MLD  and  Tj  in  Table  6.  The  mean  observed  MLD  (  H  =  60.9  m)  is 
almost  double  that  for  the  complex  (  H  =  36.5  m)  and  simple  (  H  =  36.7  m) 
events,  but  the  mean  observed  T*  (  T  s  =  8.5°C)  is  close  to  the  model-predicted 
value  for  the  complex  event  (  T  s  =  8.6°C). 


TABLE  6.  ANNUAL  MEAN  VALUES  OF  MLD,  TEMPERATURE,  NET 
SURFACE  HEAT  FLUX  AND  MIXED  LAYER  SALINITY 


Event 

Case 

Complex 

Standard 

36.5 

8.6 

-10.6 

31.9 

1 

31.9 

11.8 

-38.4 

31.9 

2 

36.8 

8.5 

-10.6 

32.0 

Simple 

Standard 

36.7 

9.7 

-21.4 

32.0 

1 

31.2 

12.5 

-43.3 

31.9 

2 

35.2 

9.7 

-21.4 

32.0 

Mean  observed  MLD  =  60.9  m  Mean  observed  Tj  =  8.5°C 


Summarizing  the  discussion  above,  we  conclude  that  the  model 
output  of  MLD  for  both  complex  and  simple  events  are  in  close  agreement 
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with  observed  MLD  in  summer.  The  estimated  mixed  layer  temperature  for 
both  events  compared  favorably  with  observations,  but  the  estimated  mixed 
layer  temperature  for  the  complex  event  is  a  little  cooler  in  summer  while 
the  estimated  mixed  layer  temperature  for  simple  event  is  higher  in  the  same 
period.  The  mean  observed  mixed  layer  temperature  is  very  close  to  that  for 
the  complex  event.  The  model  results  suggest  that  observed  data  is  best 
simulated  by  model  predicting  of  the  complex  event. 
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Figure  28.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  Tj  in  the 
Variable  Cloud  Simulation  (Simple  Event) 
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Figure  29.  Variation  of  (a)  Differential  MLD  and  (b)  Differential  Tg  in  the 
Variable  Precipitation  Simulation  in  Simple  Event 
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Figure  30.  Variation  of  (a)  Mixed  Layer  Depth  and  (b)  Mixed  Layer 
Temperature  in  Annual  Period:  A  solid  line  indicates  the  observed  data,  a 
dotted  line  indicates  the  output  of  standard  case  in  complex  event,  a  dashed 
line  indicates  the  output  of  standard  case  in  Simple  Event 
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VI.  SUMMARY 


This  study  was  conducted  to  examine  the  cloud  effect  on  the  short-term 
and  seasonal  evolution  of  ocean  mixed  layer  depth,  temperature  and  salinity 
structure  at  mid  latitudes.  The  data  for  this  study  was  collected  at  Ocean 
Station  Papa  located  in  the  Northeast  Pacific  ocean  during  1959.  Three  kinds 
of  experiments  in  association  with  different  time  scales  are  designed  for  this 
study.  Each  experiment  consists  of  three  cases.  The  Case  1  simulates  a  cloud 
forcing  with  a  20%  reduction  of  cloud  cover  from  observed  data.  The  Case  2 
simulates  a  simple  precipitation  forcing  with  E-Pr=0.  These  two  cases  are 
compared  with  a  standard  case  for  which  the  model  cloud  are  taken  as 
observed  . 

Three  different  time  scales  are  selected  for  this  study,  that  is,  3-day  period, 
30-day  and  365-day  periods.  For  a  3-day  period,  we  examine  the  albedo  and 
greenhouse  effects  of  cloud  and  the  effect  of  short-period  precipitation  events 
on  the  ocean  mixed  layer.  Simulations  using  different  amounts  of  cloud 
cover  indicate  that  the  sensitivity  of  cloud  effect  on  upper  ocean  dynamics 
depends  on  the  direction  of  the  surface  buoyancy  flux  and  the  length  of  the 
day.  In  winter,  the  upward  buoyancy  flux  is  stronger  and  the  length  of  the  day 
is  shorter,  so  the  greenhouse  effect  becomes  significant  on  the  ocean  mixed 
layer.  In  contrast,  the  downward  buoyancy  flux  is  stronger  and  the  length  of 
the  day  is  longer  during  the  summer.  Thus  the  albedo  effect  become 
significant  for  ocean  mixed  layer  dynamics.  For  simulation  using  different 
amounts  of  precipitation,  however,  we  find  that  the  ocean  mixed  laver  is 
more  sensitive  to  precipitation  in  winter  when  wind  stirring  is  predominant. 
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Similar  results  are  found  for  a  month-long  simulation.  One  striking 
feature  in  this  period  is  that  both  cloud  or  precipitation  have  more  significant 
effects  on  mixed  layer  temperature  and  salinity  when  the  mixed  layer  is  very 
shallow. 

For  long-term  “complex"  cloud  and  precipitation  simulations,  monthly 
mean  value  of  cloud  cover  and  precipitation  from  meteorological 
observations  is  derived  and  tuned  to  yield  realistic  cloud  coverage  and 
precipitation  input  for  the  model.  With  respect  to  seasonal  variation  of  solar 
radiation,  the  albedo  effect  is  significant  during  the  warming  season  when  the 
solar  radiation  is  strongest.  During  the  cooling  season,  the  longwave 
radiation  is  significant.  Therefore,  the  greenhouse  effect  plays  an  important 
role  on  upper  ocean  dynamics.  The  precipitation  effect  on  the  ocean  mixed 
layer  is  more  sensitive  during  the  cooling  season  due  to  larger  evaporation  in 
winter  and  smaller  precipitation  in  summer.  The  resulting  values  of  MLD, 
temperature  and  salinity  for  this  complex  event  are  compared  with  output 
values  for  which  the  model  is  forced  by  “simple"  constant  yearly  mean  value 
of  cloud  cover  and  precipitation.  In  general,  ine  model  output  of  MLD  in  the 
simple  event  is  shallower  than  that  in  the  complex  event.  The  mixed  layer 
temperature  in  the  simple  event  also  are  higher  than  in  the  complex  event 
during  the  warming  season.  As  for  salinity  predictions,  it  shows  that  the 
monthly  mean  values  of  observed  precipitation  and  calculated  evaporation 
(i.e.  complex  event)  are  better  used  as  initial  condition  for  the  model. 

Finally,  the  model-predicted  values  under  complex  and  simple  events  are 
compared  to  observed  BT's  values.  The  mean  observed  MLD(  H  =  60.9  m)  is 

much  greater  than  model-predicted  values  for  both  complex!  H  =  36.5  m) 
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and  simple  events(  H  =  36.7  m).  The  mean  values  of  Ts  indicates  that  the 
mean  observed  Ts(  T  s  =  8.5°C)  is  close  approximated  by  the  model-predicted 
values  for  the  complex  eventf  T  s  =  8.6°C).  It  appears  that  the  observed  MLD 
is  in  close  agreement  with  that  predicted  by  the  model  in  summer  ,  with  the 
greatest  difference  occurring  in  winter.  As  with  MLD,  the  temperature  of 
complex  and  simple  events  generally  approximate  the  temperature  cycle  and 
form  of  the  actual  data.  The  temperature  in  summer  provides  the  greatest 
difference  in  comparison  with  observations  in  the  experiment.  The  mean 
values  of  temperature  show  that  model  prediction  is  close  to  the  actual 
temperature  by  using  monthly  mean  value  of  observations. 

There  are  two  reasons  why  observed  MLD  and  temperature  may  not  agree 
with  model  predicted  value: 

1.  The  inaccuracies  in  the  BT's  data  (ship  not  at  fixed  location  &  real 
ocean  not  one-dimensional). 

2.  The  inaccuracies  in  model  forcing,  especially  in  producing  surface 
buoyancy  flux  in  variable  cloud  simulation. 

The  cloud  model  scheme  case  in  this  study  is  very  crude,  therefore 
inaccuracies  in  calculating  buoyancy  flux  for  the  model  in  variable  cloud 
experiment  m.ay  exist.  As  mentioned  before,  the  cloud  is  the  most  significant 
factor  which  influences  the  sea  surface  temperature  due  to  the  reduction  in 
solar  radiation  or  emission  of  longwave  radiation  back  to  the  ocean  surface. 
The  decreasing  or  increasing  sea  surface  temperature  may  diminish  or 
produce  more  clouds.  This  feedback  mechanism  may  couple  the  cloud-ocean 
mixed  layer  system.  In  this  study,  we  are  only  concerned  with  the 
mechanism  of  clouds  on  the  ocean  mixed  layer,  therefore  the  resulting 
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analysis  might  be  biased  by  lack  of  information  of  this  feedback  mechanism  in 
the  cloud-ocean  mixed  layer  system. 

It  is  recommended  that  this  model  experiment  should  be  expanded  by 
coupling  with  an  atmospheric  boundary  layer  model  that  includes  feedback 
mechanisms  between  atmosphere  and  ocean.  Furthermore,  the  precipitation 
induced  from  satellite  imagery  can  be  used  as  input  to  the  NFS  mixed  layer 
model  for  getting  a  more  accurate  solution.  Finally,  this  experiment  may 
help  to  focus  attention  on  the  importance  of  cloud  variability  on  ocean  mixed 
layer  because  it  has  been  shown  here  to  be  the  most  complex  and  least 
understood  aspect  of  upper  ocean  dynamics. 
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